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The thermodynamics of the classical frustrated spin chain near the transition point between 
the ferromagnetic and the helical phases is studied. The calculation of the partition and spin 

7-H ' correlation functions at low temperature limit is reduced to the quantum mechanical problem of a 

particle in potential well. It is shown that the thermodynamic quantities are universal functions of 
the temperature normalized by the chiral domain wall energy. The obtained behavior of the static 
£jjy structure factor indicates that the short-range helical-type correlations existing at low temperatures 

on the helical side of the transition point disappear at some critical temperature, defining the Lifshitz 
point. It is also shown that the low-temperature susceptibility in the helical phase near the transition 
point has a maximum at some temperature. Such behavior is in agreement with that observed in 

f"^ ' several materials described by the quantum s — 1/2 version of this model. 

m ! 

O \ I. INTRODUCTION 

^ . 

^ Strongly frustrated low-dimensional magnets attracted much attention last years A very interesting class of 

such compounds is edge-sharing chains where CuO^ plaquets are coupled by their edges 043- An important feature 

-i_> ■ of the edge-sharing chains is that the nearest-neighbor (NN) interaction J\ between Cu spins is ferromagnetic while 
the next-nearest-neighbor (NNN) interaction Ji is antiferromagnetic. The competition between them leads to the 
frustration. A minimal model describing the magnetic properties of these cuprates is so-called F-AF spin chain model 

^ the Hamiltonian of which is 

*^ . H = J\ S n ■ S„+i + J2 S„ ■ S„+2 (1) 

s : 

O . where J\ < and J 2 > 0. 

The ground state phase diagram of the s = 1/2 F-AF model has been studied extensively last years The 
ground state of this model is governed by the frustration parameter a = J 2 / \ J\\- At0<a<l/4 the ground state 
is fully ferromagnetic. At a > 1/4 the incommensurate singlet phase with short-range helical spin correlations is 
realized. At a = 1/4 the quantum phase transition occurs. Remarkably, this transition point does not depend on a 
spin value, including the classical limit s — > 00. 

However, the influence of the frustration on low temperature thermodynamics is less studied, especially in the 
vicinity of the ferromagnet-helimagnet transition point. It is of a particular importance to study this problem, 
because there are several edge-sharing cuprates with a ~ 1/4 which are of special interest [l3j . 

Unfortunately, at present the low-temperature thermodynamics of quantum s = 1/2 model Q at a ^ can be 
studied only either by using of numerical calculations of finite chains or by approximate methods. On the other hand, 
the classical version of model (JXJ) can be studied exactly at T — > 0. It is expected that the main qualitative features of 
the quantum low-temperature thermodynamics can be reproduced correctly in the framework of the classical model. 
Besides, the classical limit can be used as a starting point to study quantum effects. Therefore, the study of classical 
model |T]) is useful for the understanding of the low-temperature properties of the quantum F-AF chains. 

At zero temperature classical model fl} has magnetic long range-order (LRO) for all values of a: the ferromagnetic 
LRO at a < 1/4 and the helical-type LRO with the wave number Q = cos _1 (l/4a) for a > 1/4. It should be 
noticed that both the time-reversal and the parity symmetries are broken in the helical ordered phase. The helical 
ordered state possesses a two-fold discrete chiral degeneracy (in addition to the usual spin-rotational degeneracy) 
corresponding to clockwise and counter-clockwise turn of spins. The chiral order parameter is defined by the chirality 
vector K n = [S n x S n +i]. At finite temperature the LRO is destroyed by thermal fluctuations, the thermodynamic 
quantities have singular behavior at T — > with the corresponding critical exponents depending on a. In Ref.fhij we 
studied classical model (JTJ) at a = 1/4, i.e. exactly at the ferromagnet-helimagnet transition point. It was shown that 
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the correlation length l c behaves as T -1 / 3 and zero-field susceptibility \ diverges as T~ 4 / 3 in contrast with the ID 
Heisenberg ferromagnet (HF) (a — 0) where l c ~ T _1 and \ ~ - 

In this paper we focus on the low-temperature thermodynamics of the F-AF classical model in the vicinity of 
the ferromagnet- helimagnet transition point. As it will be shown, in the low-temperature limit the thermodynamic 
quantities become the universal functions of the scaled temperature t — T/7 3 / 2 , where the parameter 7 = (4a — 1) 
characterizes the deviation from the transition point and is assumed to be small. The calculation of the thermodynamic 
quantities in this limit reduces to the solution of the Schrodinger equation for the quantum particle in an anharmonic 
potential. 

It is worth noting that the thermodynamics of classical spin model (JTJ has been studied before [H, Oj} using the 
transfer-matrix method with the subsequent solution of the corresponding integral equations by Gaussian integration 
formulas. This method gives reliable results for the moderate temperatures and values of 7. However, this procedure 
fails for low temperatures because of appearing of an artificial gap in the excitation spectrum |18j . Our approach 
gives exact asymptotics of the thermodynamic quantities at T — > 0, 7 — > and these results are complementary to 
those obtained in [l6|, H3] • 

The paper is organized as follows. In Section II the continuum version of the model is introduced. The problem of 
the calculation of the partition function is reduced to the solution of the Schrodinger equation for a particle in the 
potential well. In Section III the spin correlation function is obtained and the fc-dependence of the static structure 
factor is studied. It is shown that the low temperature thermodynamics is the universal function of the scaling 
variable. The phase diagram of the model is constructed and the location of the Lifshitz point is determined. The 
behavior of the chiral correlation function is also studied. In Section IV the summary of the obtained results is given 
and their relation to early studies of the considered model [HI, [l?} is discussed. In the Appendix the derivation of the 
energy of chiral domain wall is given. 



II. PARTITION FUNCTION 



In the vicinity of the transition point it is convenient to rewrite Hamiltonian ([T]) in the form [7J 



(2) 



where 7 = (4a — 1). 

In Eq.([2]) we put |Ji| = 1 and omit unessential constant. In the classical approximation the spin operators S„ are 
replaced by the classical unit vectors S n . We will investigate mainly the case 7 > corresponding to the helical side of 
the ground state phase diagram. In the limit T — > the thermal fluctuations are weak and at 7 — > the period of the 
helix is long, so that neighbor spins are directed almost parallel to each other. Therefore, we can use the continuum 
approach replacing S n by the classical unit vector field s(x) with slowly varying orientations, so that 



Sn+l S n - 



ds(x n ) 
dx 



and 



S n -\ 



2S n + S n —i 



d 2 s(x n ) 
dx 2 



(3) 



(4) 



where the lattice constant is chosen as unit length. 

Using Eqs.Q and (|4| Hamiltonian ([2]) goes over into the energy functional of the vector field s(x): 



dx 



I 2 

V dx 2 



47 



d£ 2 
dx 



(5) 



This energy functional is a starting point of the investigations of model ( [TJ in the vicinity of the transition point 
a = 1/4. The partition function is a functional integral over all configurations of the vector field on a ring of length L 



"I 

It is useful to scale the spatial variable as 



Ds(x) exp < / dx 



f = 2T 1/3 a 



crs 
dx 2 
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dx 



(6) 



(7) 
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Then, the partition function takes the dimensionless form 



Z = J Ds(Z)expl-J d£ 



d 2 s\ 2 ,fds" 2 



where 
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(8) 



(9) 



T 2/3 

and the rescaled system length is A = 2T 1 / 3 L. 

Since the integration in the partition function occurs over all possible spin configurations, we are free to choose any 
local coordinate system. It is convenient to choose it so that the z axis at the point £ is directed along the spin vector 
s(£), so that the spin vector s(£) = (0,0, 1). 

Let us introduce a new vector field 

ds 

oXO = ^7 = (fe> Qv>Q*) ( 10 ) 
Then, the constraint s^(£) = 1 converts to the relations for <f(£): 

q z = 

2 _ 
iz Hx 

where the prime denotes the space derivative d/d£. 



> = -ii-ii (ii) 



Then, the first term of the Hamilton function in Eq. {[5j) transforms to 

+ # + (12) 

Here we see that the constraint s 2 — 1 effectively eliminates the q z component from the Hamilton function. 
Therefore, henceforth we deal with the q x and q y components of the vector field q only, and we denote a two-component 
vector field by q(£) ={q x ,q y )- 

The partition function in terms of q(£) takes the form: 

Z = J Dq(0 cxp | - £ d£ (q' 2 + q 4 - 6q 2 ) J (13) 

If we replace £ by an imaginary time £ — > it then partition function (I13[) takes the form of a path integral of a 
quantum particle in a potential U(q) — q 4 — 6q 2 . In other words Z in Eq.(|13p is the partition function of the quantum 
system at 'temperature' 1/A described by the Hamiltonian 

Ho = -jA + q 4 -6q 2 (14) 

where A = 9 2 + 9 2 is two-dimensional Laplace operator. Hamiltonian (fT4|) commutes with the ^-component of the 

angular momentum l z and eigenstates "0(q) of the corresponding Schrodinger equation are divided into subspaces of 
azimuthal quantum numbers l z = 0, ±1, ±2 . . . 

Ho1p n ,l, = £n,l z 1pn,l, (15) 

Thus, the wave function ip rit i z (q) describes a particle with the azimuthal quantum number l z in the 2D axially 
symmetrical potential well U(q) = q 4 — bq 2 (q = |q|). Obviously, the eigenfunctions and the eigenvalues of Eq.(fT5|) 
are real for all values of b. The dependencies of several lowest eigenvalues of Eq. (|15[) on model parameter b are shown 
in FigfT] For large b the potential U(q) takes the form of the Mexican hat with a deep and steep axially symmetrical 
well U(q) w 2b(q — a/6/2) 2 — & 2 /4. In this case the spectrum becomes almost degenerate over angular momentum l z , 
which is in accord with FigfT] In the limit b — > oo the asymptotic of eigenvalues of Eq. (|15p can be found explicitly: 

b 2 r- I 2 



n,i,=-j + «v26+^ (16) 



The partition function is a sum of exponents over all states and in the thermodynamic limit (A = 2T 1 / 3 L oo) 
only the lowest eigenvalue £o with n — l z = gives contribution to the partition function: 

z = J2e~ Xsn ^ ->e- Aso (17) 
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III. CORRELATION FUNCTION 



The calculation of the correlation function (s(l) ■ s(0)) is more complicated procedure in comparison with that for 
the partition function. It is considered in detail for the case a = 1/4 (b = 0) in Ref.[14j. where the problem of 
the calculation of the correlation functions reduces to solution of the additional differential equations together with 
Eq. (fT5|) . This procedure can be extended to the case a ^ 1/4 straightforwardly. In particular, the expressions for 
the correlation function remains exactly the same as for the case a — 1/4 and only the form of the potential well in 
the corresponding pair of differential equations is modified, U(q) = q 4 — bq 2 instead of U(q) — q 4 . The differential 
equations are 



"T"T-7 ~ ~r -J- + -r^ v + 1 v ~ bq 2 v + qu = r]v 
4 dq 2 Aq dq Aq 2 



ld 2 v 1 dv 1 
4 dq 2 Aq dq 

1 d 2 u 1 du 4 , 2 . 

—4de-T q Tq +qU ~ bqU - qV = VU (18) 

Eqs. ([18[) represent the eigenvalue problem for r} n and eigenfunctions u n (q) and v n (q) satisfying the normalization 
conditions: 

{u* n \u m ) - (Vn\v m ) = 5„, m (19) 

System of equations (fTS)) describes a two-level system in the axially symmetric potential well U(q) = q 4 bq 2 , 
where two levels with angular momenta l z — and l z = 1 are coupled by non-Hermitian transition operator. 
As shown in Ref. the spin correlation function has a form 

(s(l) ■ s(0)) = E (^\u n ) 2 e~ 2Tl/3 ^-^ (20) 

n 

where tpo and eo are the ground state wave function and the ground state energy of the Schrodinger equation ([15)) 
and u n and rj n are the eigenfunctions and the eigenvalues of Eqs. (fT8| . 

The numerical solution of Eqs. ([T8l) demonstrates the following properties. For small b all eigenvalues of Eq.([T8|) 
are real but the eigenfunctions of Eq. p^)) are complex. With the increase of b, two lowest levels 771,2 approach to 
each other, then at bo — 0.145 they coincide a nd after th at, for b > bo, these eigenvalues become complex conjugated. 
Their imaginary part increases as ^(771,2) ~ s/ (b — &o)/2. These features of the solutions of Eqs. (IT8l) are illustrated in 
Figl21 where the dependencies of the eigenvalues 771.2 on parameter b are shown. According to Eq. (|20[) the appearance 
of imaginary parts of eigenvalues 77„ implies that the correlation function begins to oscillate. 

With the further increase of b the next two levels 773,4 approach to each other and at some value b\ > bo become 
complex conjugated. And so on, so that at large values of b all low-lying levels are divided on complex conjugated 
pairs. The corresponding eigenfunctions of these pairs are complex conjugated as well. This fact guarantees the reality 
of the correlation function. 

It is important that all the eigenvalues So, r] n and the eigenfunctions ipo, u n presented in Eg. (1201) depend on the 
parameter b only. This means that the correlation function and other related physical quantities are universal functions 
of the scaling parameter b. 
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FIG. 2: Dependencies of eigenvalues 771 and 772 on parameter b. Eigenvalues 771 and 772 are real for b < 0.145: real part of 771 
(thick solid line) and rp, (thin solid line). Eigenvalues 771 and 772 are complex conjugated at b > 0.145: 3?(t7i) = 3^(772 ) (thick 
solid line) and Q^r/i) = —9(772) (dashed line). 




FIG. 3: The dependence of the normalized susceptibility X = X7 on the scaled temperature t. 



In what follows it is convenient to consider the system properties as functions of the scaling parameter t = 6~ 3 / 2 , 
which represents the normalized temperature t = T/7 3 / 2 . As shown in the Appendix the excitation energy E^w of 
a chiral domain wall, which separates two domains with different chirality vector, is proportional to E^w ~ 7 3 / 2 . 
Therefore, the parameter t is nothing but the ratio (up to numerical factor) T/Edw and the physical meaning of the 
scaling variable t is the normalization of the temperature by the energy of the chiral domain wall. 

According to Eq. ([20|) the static structure factor is 



S(k) 



— Y 



(lpo\Un) 



£0 



(Vn - £o) 2 + k 2 



(21) 



where k — fc/2T 1 / 3 is the scaled momentum. As follows from Eq. (|2"Tj) . the structure factor is scaled as T -1 / 3 and, 
therefore, S , (fc)T 1 / 3 is the universal function of t. 

The susceptibility x(k) is connected with the static structure factor by the relation x{k) — S(k)/3T. At first, we 
consider the temperature dependence of the uniform susceptibility %(0). According to Eq. (f2"Tj) it is 



^(0) 3^4/3 X] 



(lpo\Un) 2 
Vn ~ SO 



(22) 



Since the eigenfunctions and the eigenvalues of Eqs.(fl~5f and (fl"8|) are functions of a single variable t, the normalized 
susceptibility x = X7 2 can be expressed as an universal function of the normalized temperature t. Thus, x reveals 
the scaling behavior: in the vicinity of the transition point and at T — > the dependence of x on the temperature 
and the model parameter 7 enters only in the combination T/7 3 / 2 . The dependence of x on t is shown in Figj3] The 
characteristic features of this dependence are the existence of maximum of x at t = t m (t m ~ 0.432) and the finite 
value of x at T = 0. The latter fact is a pure classical effect. In the s = 1/2 quantum F-AF model the ground state is 
believed to be gapped (though the gap can be extremely small [7|), and so x -> as T -> 0. The obtained dependence 
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FIG. 4: The static structure factor S(k) for several values of the scaled temperature t. 



x(t) means that for a given model parameter 7 the uniform susceptibility has the maximum at temperature T m ~ 7 s / 2 
and the height of this maximum is x-m ~ 7 -2 - This implies that with the increase of 7 the maximum of x shifts to 
higher temperatures and the magnitude of the maximum x-m decreases. 

Let us consider the behavior of the static structure factor (|2 1 [) as a function of k. Some examples of S(k) are shown 
in Fig|31 As we see, the static structure factor has a pronounced maximum and its location depends on t. In the 
vicinity of the maximum S(k) can be represented as 



S(k) = 



S(k m ) 



i + p c (k-k m y 



(23) 



Therefore, the behavior of the static structure factor is described by three parameters: the location of the maximum 
k m , the magnitude of the maximum S(k rn ), and the correlation length l c characterizing the width of the maximum: 



l.=- 



1 



(24) 



^S\kjff\ dk^ 

Obviously, the correlation length is scaled as the space variable and, therefore, the normalized correlation length 
l c = 2T 1 / 3 ^ C is the universal function of t. 

The location of the maximum k m is determined by the condition dS/dk = 0. According to Eq. (|21[) this condition 
reduces to the equation 



~ k (V'ol"") 2 (Vn - £0) = Q 



(25) 



This equation determines the scaled momentum k m as a function of t. Obviously, Eq. (|25l) always has a trivial 
solution k m = 0. For high values of t there are no other solutions. However, at some value of the scaled temperature 
t = t c another, non-trivial solution of Eq. (125[) with k m ^ appears. For t < t c the maximum of the structure factor 
shifts from k m = 0. Thus, t — t c defines the Lifshitz point which corresponds to the transition from the helical to 
the ferromagnetic phase. Certainly, there is no true helical or ferromagnetic LRO in these phases, because thermal 
fluctuations destroy any LRO in one-dimensional systems. Therefore, under the helical and the ferromagnetic phases 
we mean the presence of the short-range order (SRO) of the corresponding type. 

The location of the Lifshitz point is determined by the equation: 



E 



(Vn - £0) 



(26) 



The numerical solution of Eq. (|26|) gives t c ~ 0.925. Thus, the transition line between the ferromagnetic and the 
helical phases ('Lifshitz boundary') in (T, 7) plane has a form 



T c ( 7 ) = 0.925 7 



3/2 



(27) 



The corresponding phase diagram is shown in FigJS] 
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FIG. 5: The phase diagram of model (0). 
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FIG. 6: The dependence of k m /Q (Q = \f2"f) on the scaled temperature t. 

For a fixed model parameter 7 < 1 the dependence of the system properties on the temperature can be qualitatively 
described in the following way. At T = the structure factor S(k) has a (5-function form at the corresponding wave 
vectors k m — ±t/2j, indicating the helical LRO. At low temperatures (T <C T c ) the 5-peaks are smeared by gapless 
excitations (spin waves) over the helical ground states. The chiral domain wall excitations have the gap ~ 7 3 ^ 2 ) 
and their contribution to the thermodynamics is exponentially small. The analysis of Eq. (|21[) shows that in this region 
(T <C T c ) the magnitude of the maximum and the correlation length behave as S(k rn ) ~ l c ~ "f/T. 

When the temperature is increased, the maximum of S(k) shifts to k = and at approaching to the Lifshitz point 
the value of k m tends to zero by the square root law k m ~ r j~ 1 ^ i \/T c — T. This is similar to the order parameter 
behavior in Landau theory of second-order phase transitions. Indeed, in the three-dimensional frustrated classical 
spin model, when the helical LRO exists, the wave vector associated with the helical phase vanishes by the same law 

5 T—1 1 
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FIG. 7: The dependence of the normalized correlation length l c — 2T 1//3 Z C on the scaled temperature t. 
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20]. 

The numerically obtained dependencies of k m /^2rj and the scaled correlation length l c on t are demonstrated in 
Fig|6] and FigJTl correspondingly. As shown in Fig[7] the correlation length vanishes at t = t c , because S(k) does 
not contain k 2 term in k expansion at the Lifshitz point. The magnitude of the maximum at the Lifshitz point is 
5(0) « I.97- 1 / 2 . 

For T » T c the helical correlations are destroyed, the maximum of S(k) is situated at k rn — and in the limit 
t -> 00 it tends to the results obtained in Ref.pJ]: S(0) = 3.21/T 1 / 3 and l c = 1.04/T 1 / 3 . 

Up to now we considered the classical spin model where S is usual three-component classical vector. It is appropriate 
to discuss briefly the planar model describing by Hamiltonian (|5J), where the vector S has two components (planar 
spin model). For short we denote below these two models as v = 3 and v = 2 cases. Just as for the v = 3 model the 
transition point between the ferromagnetic and the helical ground states for the v = 2 model takes place at a = 1/4. 
The low-temperature thermodynamics of the v = 2 model near the transition point can be analyzed in a similar 
manner as for the v = 3 case. In particular, the calculation of the partition function for the v = 2 model reduces to 
the solution of the Schrodinger equation 



■71%+ U b)4> = e1> (28) 



where U(q) = q A — bq 2 . 

Our study of the planar model shows that the thermodynamic quantities are the universal functions of the same 
scaling parameter t — T/7 3 / 2 . The behaviors of S{k), k m (t) and l c (t) are similar to that for v = 3 case. The Lifshitz 
point is determined by t c ~ 0.7 and the critical exponent characterizing the behavior of k m near t c is 1/2 as well. 
However, the behavior of the chiral correlation function 



R l = {K n -K n+ i) (29) 
of these two models are different. In the continuum approach the chiral vector is expressed as 

K n = S n x S n+1 w s(x) x ^ (30) 
In the v = 2 model the chiral vector is directed perpendicular to the spin plane for any x. Therefore in this case 

^)^(o)«^-« = #Mo) (31) 

In the v = 3 model we take into account local coordinates with s = (0, 0, 1) and q z — and obtain similar expression 

K(l)-K(0)=d(l)-q(0) (32) 

The expression for the chiral correlation function Ri can be obtained in the same manner as for the spin correlation 
function: 

^=£<^o|q|^} 2 e - 2Tl/3(£ "- £o) ' (33) 
n 

Here ip n and e n are the eigcnfunctions and the eigenvalues of the Schrodinger equations (fT5|) . For planar case Ri has 
similar expression where ip n and e n are solutions of Eq. (|28l) . Since ip n and e n are always real, the chiral correlation 
function Ri does not oscillate and the maximum of the corresponding structure factor is at k — for both v — 2 and 
v = 3 cases. 

As follows from Eq. (|3"3"|) the correlation length of the chirality £ is scaled as T 1 ' 3 and at T — > it is defined by the 
lowest level having the non-zero matrix element with the ground-state wave function ipQ. So, for normalized chiral 
correlation length ( = 2T 1 / 3 £ we have 

C=-^- (34) 

£l — Eq 

The ground state wave function for Eq. (|28[) is an even function of q and only the eigenfunctions with odd parity 
give contributions to i?/. So, for the v = 2 model £1 is the energy of the first excited state of the odd parity. For the 
v = 3 model the matrix elements in Eg. (1331) are non-zero only for the wave functions ip n with l z = 1. At the transition 



9 



point (corresponding to the limit t — >• oo) the energy difference Si — £o is finite for both v — 2 (e% — Eq — 1.09) and 
v = 3 (ex — £q = 1-21) cases. 

At t — >• the behavior of the energy gap between two lowest levels for v = 2 and v = 3 models becomes qualitatively 
different. For the v = 2 case the potential £/ (g) = q 4 — bq 2 has two degenerate minima and the splitting of two lowest 
levels £\— So depends upon tunneling through the classically forbidden region of small q. This tunneling is exponentially 
small for low temperatures and the WKB estimate gives £ ~ exp(l/3-\/2*)- However, for the v = 3 case the particle 
moves in two-dimensional axially symmetric potential, the form of the Mexican hat, and relevant splitting is equal 
to the additional energy associated with one unit of the angular momentum. Then, according to Eq. flTSl) C — 4/t 2 / 3 . 
Thus, the behavior of the correlation length of the chirality is essentially different for these two models. 

So far we studied the helical side of the transition point a > 1/4 (7 > 0). On the ferromagnetic side (7 < 0) the 
static structure factor S(k) and the susceptibility x(^) have the maximum at k — for any temperature. The uniform 
susceptibility x diverges at T — > but the corresponding exponent changes from 4/3 to 2 and Eq.(f!?2"j) describes 
this crossover. The exponent 2 appears in the limit b — > —00, when = — 26 3 / 3 and the susceptibility becomes 
X = 2|7|/3T 2 , which is in accord with the result for HF with the coupling constant |7|. Indeed, in the limit b — > —00 
one can neglect the second derivative term in energy functional ([SJ and the model reduces to the HF with renormalized 
coupling |7|. 

IV. SUMMARY 

We have studied the low-temperature thermodynamics of the classical F- AF model in the vicinity of the ferromagnet- 
helimagnet transition point with use of the continuum approximation. The calculation of the partition and the spin 
correlation functions is reduced to a quantum mechanical problem of a particle in a potential well. It is shown that 
the thermodynamic quantities are universal functions of the scaling variable t = T/7 3 / 2 , where 7 = 4J2/IJ1I — 1 
describes the deviation from the transition point. 

On the ferromagnetic side (7 < 0) the static structure factor S(k) has the maximum at k = for all temperatures. 
The uniform susceptibility \ diverges at T — > but the corresponding asymptotics changes from x ~ T~ 4 / 3 at the 
transition point to \ ~ \^l\/T 2 for T -C |7| 3 ^ 2 - The crossover between these regimes takes place at T ~ 7 s / 2 . 

On the helical side of the transition point (0 < 7 < 1) our investigation displays the following dependence of the 
system properties on the temperature. At zero temperature the static structure factor S(k) has a 5-function form at 
the corresponding wave vectors k m — ±y^y, indicating the LRO of the helical-type. At finite temperatures the LRO is 
destroyed by thermal fluctuations, which is manifested in smearing of <5-peaks of the static structure factor by gapless 
excitations (spin waves) over the helical ground states. The chiral domain wall excitations have the gap (-Edw ~ 7 3 ^ 2 ) 
and their contribution to the thermodynamics is exponentially small for low temperatures T <C 7 3 ^ 2 . In this region 
the magnitude of the maximum of the structure factor and the correlation length behave as S(k m ) ~ l c ~ 7/J 1 . 

Further increase of the temperature is accompanied by the damping of the peak and its shift to k = 0. At 
approaching to the Lifshitz point T — >• T c = 0.9257 3 / 2 the value of k m tends to zero by the square root law k m ~ 
7~ 1 / 4 \/X'c — T. Such dependence is similar to the order parameter behavior in the Landau theory of the second-order 
phase transitions. So, the location of the peak k m plays a role of the order parameter, though one should remember 
that the system has only the short-range helical order at finite temperatures. Equation T c — 0.9257 3 / 2 determines 
the Lifshitz boundary in (T, 7) plane, separating the ferromagnetic and the helical phases. 

For T > T c the helical correlations are destroyed by the chiral domain wall excitations and the maximum of the 
static structure factor S(k) is situated at k rn = 0. The high temperature limit T 3> T c is equivalent to the limit 
7 <C T 2 / 3 . This implies that one can neglect the second terms in Eqs.([2]) and §5§ and the model effectively reduces 
to that at the transition point studied in Ref.[l4j]. Here the uniform susceptibility and the correlation length are 
X(0) = 1.07/T 4 / 3 and l c — 1.04/T 1 / 3 . We should remind that our analysis is restricted by low temperatures and, 
therefore, the latter high-T limit is in fact restricted by T c <C T <C 1. 

The uniform susceptibility \ has the maximum at T m ~ 7 3 / 2 and Xm ~ 7 -2 > i-e- with the increase of the frustration 
parameter a the maximum of x shifts to higher temperatures and the magnitude of the maximum Xm decreases. The 
obtained dependence of v on T is in a qualitative agreement with those observed in the edge-shared compounds with a 
close to 1/4 (LiiCuOi 2], RbiCuiMo^O\i Q, Li 2 CuZr04 [III). The dependencies of T m and Xm on the frustration 
parameter a are also in accord with the experimental observations (l3| . 

We have also shown that the thermodynamic behavior of both planar and classical F-AF models are similar. The 
exception is the chiral correlation function. The chiral correlation length of the planar model increases exponentially 
at T in contrast with the power dependence for the classical v = 3 model. This fact is explained by different 
topology of the potentials in the corresponding Schrodinger equations. 

It is interesting to compare our results with those obtained in Refs.0, [l3] for intermediate values of T and 7. In 
contrast with the phase diagram obtained for small values of T and 7 and schematically shown in Fig. 6, the Lifshitz 
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boundary for intermediate values of T and 7 can have the reentrant behavior [l6[ . There is a critical value of 7 = j c 
(for the planar model 7 C ~ 0.16 [16j). where the Lifshitz boundary turns round and moves back to 7 = with the 
further increase of temperature. This means that for a fixed small value of 7 (0 < 7 < 7 C ) the ferromagnetic phase 
exists in the finite range of temperatures and for 7 > j c the helical-type correlations remain for all temperatures. 
This fact testifies that in the region of intermediate or large values of T and 7, where the continuum approximation 
becomes unapplicable, the thermodynamic behavior can be qualitatively different from the results presented here for 
small T and 7. 



Appendix 

In the appendix we calculate the chiral domain wall energy for the classical spin model described by energy functional 
([5]). To study excitations in this model it is useful to make the following rescaling of the space variable x = £/y/2rf. 
Then, energy functional ([5]) takes the form: 



E = 



7 3/2 
23/2 



r/£ 



d^ 2 J \d£ 



(A.l) 



Here we see that the integrand in Eq. (|A.ip does not contain any parameter and, therefore, any localized excitation 



is proportional to 



r 



3/2 



Certainly, the determination of the numerical factor before 7 s / 2 needs the solution of 



the corresponding Euler equation. These Euler equations are highly non-linear and its complete analysis is hardly 
possible. Fortunately, such an analysis is possible for the planar spin case. The solutions for the planar spin case 
represent particular cases of solutions for the original three-component spin model. 

For the planar model the two-component spin vector can be represented as s = (sin 9, cos 9) and energy functional 
(|A.1[) in terms of 9(x) has a form 



E = 



7 



3/2 



2 3 / 2 



cP9_ 



- 2 



Variation of energy functional (|A.2|) in leads to the Euler equation 



d 4 9 



6 



d 2 



89 



t &8 







(A.2) 



(A.3) 



Here we see that the Euler equation has a trivial helical solution = a(£ — £0) with any constants a and £o- 
However, the minimum of energy (IA.2I) corresponding to the ground state is achieved by a definite helical configuration 
with 



9(0 = ±(£ - £ ) 



(A.4) 



and the ground state energy is 



En = —L-!— 



(A.5) 



To find the excitations over the helical ground state we should look for other solutions of Eq. (|A.3[) . Eq. (|A.3|l can 
be integrated so that for = 89/ d£, we obtain equation 



y' 2 = y 4 -2y 2 



Ay + B 



with some integration constants A and B. 

The energy functional in terms of takes the form: 



E-En = 



r 



3/2 



2 3 /2 



|d£ [y' 2 + (l-y 2 ) 2 ] 



(A.6) 



(A.7) 



The right-hand side of Eq. (|A.6[> has a double- well form. The complete analysis of solutions of Eq. (jA.6[> will be 
given elsewhere. Here we present the solution describing the chiral domain wall only. This case corresponds to the 
choice of integration constants ^4 = and B = 1. In this case the solution of Eq. (|A.6[) has a simple form 



y(£) = tanh (£ - £„) 



(A. 
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or in original terms for angle 9{x) 

6(x) = d + In cosh U/2y(x - x J (A.9) 

with arbitrary 9q and xq- As can be seen, this solution describes a chiral domain wall separating two domains of 
opposite chirality. The energy of this excitation is 

£dw = ^(2 7 ) 3/2 (A.10) 

In case of three-component spin vectors the chiral domain wall excitations are more complicated. They depend on 
the boundary conditions on x — ±00, which can be imposed on the chiral vectors K. In the helical ground state the 
length of the chiral vectors equals \K\ = y/^rj, but the direction is arbitrary according to degeneracy of the helical 
ground state. Therefore, the chiral domain wall solution is a function of the angle between the chiral vectors at 
x = ±00: 

/ \K-qo ■ Koo) I . -1 -1 \ 

cos0= (A. 11) 

2 7 

In the planar spin case considered above the angle (j) is <f> = 7r, because chiral vectors are antiparallel on ±00. 
It is also obvious that the limit <f> — > corresponds to the ground state configuration. In general the dependence 
of the energy of the chiral domain wall on the angle <f> is unknown, but it is natural to assume that this function 
monotonically increases from zero at <j> = to the exact result Eq. (|A.10|) at <\> = it. Therefore, in any case the energy 
of the chiral domain wall is proportional to ~ -y 3 / 2 and can be written as 

£dw = f(4>h 3/2 (A.12) 

with some smooth monotonic function /(</>) having the limits /(0) = and f(n) = 2 3 / 2 /3. 
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